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We study heavy-quarkonia suppression under final-state multiple scatterings in 
most central Au + Au collisions at RHIC energy. We first calculate the survival 
probability of a heavy quarkonium under multiple scattering in Bjorken's expanding 
QGP at large N c . Then, we calculate the rapidity dependence of the nuclear modifi- 
cation factor Raa for heavy-quarkonia production by considering final-state multiple 
scatterings in most central Au + Au collisions in a simplified model. In our formula a 
constant Po is also introduced to estimate the possible cold nuclear effects. By fitting 
the data for J/\I/ production in most central Au + Au collisions at ^snn = 200 GeV 
at RHIC, we find that the transportation coefficient qo — (0.33 — 0.95) GeV 2 /fm, 
and, accordingly, the energy density at tq is eo — (1.39 — 5.62) GeV/fm 3 in perturba- 
tive thermal QCD. A better understanding of cold nuclear effects is essential for us 
to get a more accurate analysis. The small values of the transportation coefficient go 
in our estimate are in sharp contrast with those obtained by the analysis of high-p^ 
hadron spectra in Ref. [31]. 
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I. INTRODUCTION 

Relativistic heavy-ion collisions at RHIC and the forthcoming LHC provide the oppor- 
tunity to study QCD matter of extreme temperature and density, the quark-gluon plasma 
(QGP)[l|. Light quarks and gluons are expected to be thermalized to form QGP in a short 
time r ~ 1 fm. Heavy quarks are too heavy to be thermalized in such a short time. With a 
formation time much shorter than To they nearly live through all the stages of heavy ion col- 
lisions and, in this way, act as an excellent probe to study the properties of the background 
medium. Because of their small size, heavy quarkonia, the bound states of heavy quarks, 
play an important role in understanding the formation and properties of the QCD matter 
created in heavy ion collisions 

J/\l/ suppression was originally proposed by Matsui and Satz as a signature for the QGP 
formation in heavy ion collisions [3|]: in QGP color screening dissolves J/\l/ into open charmed 
mesons D and D. The NA38 collaboration at the CERN SPS observed for the first time 
the suppression of J/^ production in oxygen-uranium reactions at 200 GeV/nucleon 4j. 
However, the answer to whether J/\l/ suppression can be taken as an unambiguous signature 
for the QGP formation or not becomes complicated since one can not rule out alternative 
explanations by cold nuclear effects (see and references therein). 

J/\l/ suppression has been observed with puzzles in Au + Au collisions at ^s^n — 
200 GeV by the PHENIX experiment at RHIcfl. One puzzle is that the suppression at 
midrapidty is similar to that observed at the SPS and smaller than expectations based 
on the analysis of the local medium density. This possibly implies that the enhancement 
of J/^f production by recombination of initially uncorrelated cc paires are important at 
RHIC energy (3, Isj]. The other puzzle lies in the dependence of the px integrated nuclear 
modification factor Raa on Np aT ^, the number of participant. Raa is more suppressed at 
forward rapidity than at midrapidity. Possible solutions for such a puzzle have been proposed 
in (a-ll|. 



In this paper, we study the dissociation of heavy quarkonia by final-state multiple scat- 
terings in Bjorken's expanding QGP. At RHIC energy, the heavy quarkonia produced in 
heavy ion collisions are typically non-relativisitc. Therefore, one can estimate the creation 
time Tb by the inverse of the binding energy Eb in the quarkonium rest frame. If the back- 
ground medium is thermalized within the thermalization time tq ~ 1 fm, heavy quarkonia 
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will be created before thermalization. Therefore, one can use heavy quarkonia to probe the 
thermodynamic properties of the system at r since the information about thermalization 
can survive final-state multiple scatterings. 

T, the lowest bb bound state, is a more ideal probe to decipher information about ther- 
malization of the background medium than J/ty. T, with a size (ob) — 0.28 fm, is the most 
tightly bound state and the background QCD matter, thermalized or not, is more partonic 
to such a hard probe. Due to their short formation time r# ~ 0.18 fm, T can actually 
witness and record the thermalization of the partonic background medium. Moreover, at 

n 

RHIC energy the recombination effect is small in T production |8|] such that the information 
of thermalization carried by T can survive final-state interactions. On the experimental 
side, the study of T production in relativistic heavy ion collision for the first time becomes 



possible 



121 ]. Because of the increased cross-seciont of T production at LHC energy, one can 



expect extensive study of T suppression in heavy ion collisions at LHC. 

In this paper we focus on the significance of final-state multiple scattering to heavy 
quarkonia suppression and neglect other possibly important effects, such as partonic ener 
loss 13 1 . the gluon saturation effects j^], the nuclear absorption and recombination effect[7|, 



14|. This paper is organized as follows. In Sec. [TTJ we first give a brief review about 



Bjorken's hydrodynamic model and a new derivation of the multiple scattering of a quark in 
it. Then, we calculate the survival probability of a heavy quarkonium in such an expanding 
QGP at large N c . In Sec. IHIl we investigate the implications of our formula of the survival 
probability to heavy-quarkonia suppression at RHIC in a simplified model. First, we give an 
intuitive explanation to the experimental results about J/ty Raa versus y in most central 
Au + Au collisions at y/s^N = 200 GeV. Then, we illustrate in the meanwhile how one 
can possibly abstract %, the transportation coefficient at tq by studying J/\I/ suppression in 
Bjorken's expanding QGP. This quantity is of great significance in our understanding of jet 



quenching in heavy-ion collisions 



131 ] . Conclusions are presented in Sec. [TV] 



II. DISSOCIATION OF HEAVY QUARKONIA BY MULTIPLE SCATTERINGS 

IN AN EXPANDING QGP 

In this section, we give a brief discussion about heavy-quarkonia dissociation due to 
final-state multiple scatterings in an expanding QGP. As illustrated in Fig. [1] a heavy 
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FIG. 1: Multiple scattering of a heavy quarkonium. 
quarkonium with a size (ob), after creation, travels through a slab of expanding QGP and 



we wil 



calculate the survival probability for such a process. We use Bjorken's hydrodynamic 



model 15] to describe the evolution of the background medium after the thermalization time 

T - 



A. Bjorken's hydrodynamic model 15] 



Imagine that after two gold nuclei collide head-on, the system reaches local thermody- 
namic equilibrium at time To- Following Bjorken, we assume cylindrical symmetry in the 
transverse direction to the collision axis. The energy density e and the entropy density s only 
depend on the invariant time r = \ft 2 — z 2 and respectively satisfy the following equations 



de 
dr~ 



P 



T 



and 



(1) 
(2) 



ds s 
dr t ' 

where P is the pressure. In the following, we assume that the pressure and the energy 
density are related by the following equation 



P = c% 



(3) 



with the sound velocity c s a constant. Inserting into Equ. ([T]) the following thermodynamic 
relations 



-r— = s and e = Ts — P, 
dT 



(4) 
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FIG. 2: One scattering of a quark off the QGP constituents: X is the relative transverse coordinates 
of the two quark lines. 



one obtains 



— - -r 2 - ft) 



Accordingly, one can easily get the so-called Bjorken's scaling solution 



T = T pV\ (6) 



T 



T 

and 



(8) 



T J 

where so, eo and r are respectively the entropy density, the energy density and the temper- 
ature at the thermalization time tq. 

The conservation of entropy (see Equ. ©) predicts the central-plateau structure of 

n 

particle multiplicity in the rapidity spectra, which is not observed at RHIC[16[. However, 
the deviations from boost invariance are not very large for \y\ < 2.5. Moreover, our picture 
about heavy-quarkonia suppression in the following sections can be easily generalized to 
other more complicated hydrodynamic models. 



B. The multiple scattering of a quark in an expanding QGP 17|] 



When a quark travels in the expanding QGP cylinder, it picks up a transverse momentum 
from random kicks by the QGP constituents. In the following, we will give a new deriva- 
tion of the evolution of -$^j-(p±,t,tq), the transverse momentum distribution of the quark, 
from initial invariant time r up to time r. The fundamental assumption here is that the 
quark undergoes uncorrelated multiple scatterings with the QGP constituents. Let us define 
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/(X, r, r ) such that 



dN 
d 2 p± 



(p±,r,r ) = J d 2 xe- ip - x f(X,r, 



To 



(9) 



Given /(X, t, To), one can choose an infinitesimal time interval dr such that the difference 
between /(X, t + dr, t ) and /(X, t, t ) can be calculated by considering only one extra 
scattering. Including all the three Feynman diagrams in Fig. [21 in eikonal approximation 
we have 



/(X, t + dr, t ) - /(X, t, t ) = /(X, t, T )drp(r) J 
or, equivalently, 



d 2 q da 
(2tt) 2 ^ 



(g 2 ,T)(e^ x -l), (10) 



J^(X,t,t ) = p(t)/(X,t,t ) J d 2 q^(q 2 ,r)(e^-l), (11) 



where ^(g 2 , r) = a s C F f A\V(q 2 , r)\\ 



q u 7 / -o - - j (2tt) 2 

Solving ffTTT) with the initial condition /(X, t ), we obtain 



/(X,T,T ) = /(X,T )exp|£rfT'p(T') J ^J^V)(e' q X -l)}. (12) 



The transport coefficient $(t) is defined as 18| 



q( T ) = p( r ) / d 2 qq 2 



da . 2 
<i 2 g 



For X 2 ~ 0, we have [171] 



/(X.r.m) :/(X.7- l) )ox 1 )^-i<rf.(r,7 il )>Jf 2 



(13) 



(14) 



where the pr-broadening of the quark is given by 

(A(r,r ))= f dr'q(r'). (15) 
The leading term of the high-temperature expansion for q(r) takes the following form [it| 

q(r) = q (^y 4 , (16) 

and, accordingly, 

(pl(r,r )) = q(r)r ^ ^ . (17) 
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FIG. 3: One scattering of a color singlet dipole off the QGP constituents: xx and x'j_ are respec- 
tively the relative transverse coordinates of the quark and antiquark in the amplitude and in the 
conjugate amplitude. 

C. Survival probability of a heavy quarkonium in an expanding QGP 



After creation, a heavy quarkonium travels through a slab of QGP. If neglecting higher 
Fock components, one can take the heavy quarkonium as a QQ color singlet dipole. A 
detailed discussion of the initial-state nuclear effect and the validity of the dipole model in 
describing heavy-quarkonia production in heavy- ion collisions is presented in Ref. [9J. In this 
paper, we only deal with the dissociation of a color-singlet dipole(quarkonium) after creation, 
which differs with Ref. As illustrated in Fig. [U the color singlet dipole may break up 
due to final-state interactions in the QGP. If r, the time for the quarkonium traveling in 
the QGP, is smaller than the formation time Tf orm = 'yTB, one can neglect the interaction 
between the dipole and only consider the multiple scattering. However, if r > tf orm , one 
has to consider the mutual interaction between the quark and anti-quark in the quarkonium 
besides the p^-broadening of them under multiple scattering. In this subsection, we will 
calculate the survival probability in these two cases. 



1. The multiple scattering of a heavy- quarkonium with t < t 



form 



Within the formation time 775 = 7/ 1 Eb with Eb the binding energy, one can neglect the 
interaction between the quark and anti-qaurk in a bound state [l9f] . First, let us calculate the 
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r— evolution of f(x±, x' ± , t, t ), the product of the quarkonium state in the amplitude and 
that in the conjugate amplitude. As shown in Fig. [31 xi and x^_ are respectively the relative 
transverse coordinates of the quark and antiquark in the amplitude and in the conjugate 
amplitude. In eikonal approximation, only the relative transverse coordinates of the quark 
and antiquark are relevant 20] . In the following we will take the light-cone wave function of 
the quarkonium as a function only of relative transverse coordinates. Before the initial time 
t , when there is no QGP present, we have 

f(x±, x' ± , t ) = (p(x±)(p*(x' x ), (18) 

with (f(x±) the light-cone wave function of the quarkonium. Neglecting terms proportional 
to 0(1/N 2 ), for one scattering of the heavy quarkonium off the QGP constituents we need 
only sum over the six diagrams as shown in Fig. (3^20]. Following the derivation of /(X, r, tq) 
for the multiple scattering of a quark in the previous subsection, it is easy to show that in 
eikonal approximation f(x±, x' ± , r, r ) satisfies the following partial differential equation 

-^-f(x ± ,x' ± ,T,r ) = p(r)f(x x ,x' x ,T,T ) J d 2 q^q 2 ,r){e^ + e-^-2), (19) 

where jg- q (q 2 ,r) = a s C F j j0f\V{q 2 , r)| 2 . The solution of Equ. flU]) with the initial condi- 
tion ( {TBI takes the following form 

f(x ± , x' ± , r, r ) = ¥>(*.i>V-l) exp f drp(r) J d 2 q^(q 2 , r) (e^ + e~*»< - 2 ) 

~ ^(x ± )^(x / ± )e-3<pi(-^o))(4+-'x 2 ) ; (20) 

where (pj_( T , T o)) is given by Equ. ( IT71) and we have made the approximation that x± and 
x' ± are much less than the inverse of the typical momentum transfer. 

Now, we are ready to calculate the survival probability of the heavy quarkonium up to 
time r, which is given by 

P SU r(r, r ) = J d 2 x±d 2 x' ± f(x±, x' ± , t, T )f*(x±, x' ± , T ). (21) 

According to (1201) . we have 

Psur(r,r ) = [ dV^k(^)|>(^)| 2 e-^-^^ + ^. (22) 



Let us approximate the heavy-quarkonium overlap function with Gaussian wave packet 
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TABLE I: Static properties of heavy quarkonia: in the rest frame of the quarkonium, we take the 
creation time tb = 1/Eb with Eb the binding energy and the typical momentum p = l/(a#) with 
(as) the size of the quarkonium. 



quarkonia 


mass (GeV) 


E B (GeV) 


t b (fm) 


p (GeV) 


a B (fm) 


J/¥ 


3.10 


0.64 


0.31 


0.40 


0.5 


T 


9.46 


1.10 


0.18 


0.71 


0.28 



where (as) 2 = \a 2 B with a B the size of the quarkonium, is listed in Table. H|2l|. Inserting 
Equ. fl23l) into Equ. fl22|) . after some algebra we obtain 



P S ur{T, T ) 



16 



(24) 



{A+(a B ) 2 (pl(r,T )}} 2 - 
Equ. fl22|) justifies the following physical picture for the dissociation of a heavy quarko- 
nium due to multiple scatterings 191]: if the quark and anti-quark in the quarkonium pick 
up transverse momenta greater than the typical momentum in the bound state, that is, 
a/ (p 2 ± ) > l/(a B ), the quarkonium breaks up. Since the dominant contribution to P sur 
comes from the integration of x± and x' ± over the region of x\ < (as) 2 and x' 2 < (as) 2 , one 
can easily see that 

• if (p 2 ± ) > l/{a B } 2 , P sur < 1; 

• if (pi) « l/(a B ) 2 , P sur ~ 1. 

Such a quantitative behavior analysis tells us that the transition between large and small 
survival probability must happen at (pj_) ~ 1/ (cib) 2 , and this is exactly the manifest of the 
above physical picture. 



2. The multiple scattering of a heavy- quarkonium with t > tf orm 

In this subsection, we calculate the survival probability of a quarkonium traveling in a 
slab of expanding QGP within a time t = L > tf orm . In this case, one can not neglect 
the interaction(attraction) between the quark and anit-quark. The wave function of the 
quarkonium, in the interaction picture, completely changes and the x± coordinate is not 



10 



t = 



V Q 



-8- 

o 
o 
o 
o 
o 
c 



o 
o 
o 
o 
o 



o 

~o - 

o 
o 
o 
o 
o 
c 



2 YT B 



o 
o 
o 
o 
o 
c 



o 

~o _ 

o 
o 
o 
o 
o 
c 



3YT B 



o 
o 
o 
o 
o 
o 



o 
o 
o 
o 
o 
o 



o 
o 
o 
o 
o 
c 



FIG. 4: The picture of the dissociation of a heavy-quarkonium under multiple scattering within a 
time interval t = L > tf orm = ^tb = 7 /Eg: the quark(anit-quark) is pulled once by the color field 
of anti-quark(quark) within a period At < 7T#. Since the quarkonium can already dissociate with 
some probability within At ~ tf orm , as calculated in Equ. (|22p . the mutual attraction between 
the quark and anti-quark is not enough for the binding. And the final-state survival probability 
is the product of that in each time interval of At ~ tf orm within the time interval t under the 
assumption that the only effect of interaction between the quark and anti-quark within At ~ jtb 
is to completely repair the quarkonium.. 

"frozen" any more. The complete calculation is beyond the eikonal approximation and we 
use an approximate picture illustrated in FigJUto estimate the multiple scattering effect. 

In this picture, the quark (anit- quark) is pulled once by the color field of anti-quark(quark) 
within a period At ~ 775 = j/Eb- In the vacuum, the quark ( ant i- quark) can always be 
pulled back into the bound state, that is, the quarkonium can completely repair itself. In 
the medium, because of screened inter-quark potential and the multiple screening effect, the 
mutual attraction between the quark and anti-quark might not be enough for the binding. 
Still, we assume that the only effect of interaction between the quark and anti-quark within 
At ~ 7773 is to completely repair the quarkonium. Obviously, under such an assumption we 
totally neglect the color screening effect and we will give a more detailed discussion in the 
Conclusion. The picture could give the parametrically right answer b y ta king At ~ 7T# = 



'j/Eb even though there is always a constant we can not really control 19j. In the following, 
we take At = £775 with £ a constant. Therefore, the final-state survival probability is the 
product of that in each time interval of £775 within the time interval t, that is, 

(i~,r n ), (25) 

where the time for the quarkonium travels in the medium is t < (n + 1)^tb and Tj is the 
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invariant time corresponding to i^Ts- 
In perturbative thermal QCD, one has 



c 2 s = \ + 0{al). (26) 



Keeping only leading terms in 0(a s ), we have 



T 



<P±(t,t )) ~q r \n-, (27) 
To 

Taking the QGP as aweekly coupled gluon gas, q and the energy density e respectively take 
the following forms 



22 



q = ^.a^T 3 , (28) 

7T 

and 

e = -j^-T 4 , (29) 

with C(3) ~ 1.202. 

In the following sections, we will only take Tq and %{To) as parameters and choose (? s = |. 
The transportation coefficient go characterizes the thermodynamic properties of the back- 
ground medium at the thermalization time tq. With tq and go fixed by RHIC's data about 
J/ty suppression, it is possible for us to decipher some information about thermalization. 



III. J/* SUPPRESSION AT RHIC 

In this section we investigate the significance of final-state multiple scatterings to heavy- 
quarkonia suppression in a simplified model illustrated in Fig. [5j We assume that the 
background medium is thermalized in a short time r ~ 1 fm < Ra, the radius of the 
QGP cylinder. From Table [1^21 ] . one can see that the formation time of heavy quarkonia 



tb is shorter than r . To simplify our discussion, we use the Glauber model [23j at the stage 
of heavy-quarkonia creation, that is, the heavy quarkonia are created in nucleon-nucleon 
binary collisions of the two colliding nuclei at time tb- After tq, the quarkonia travel in the 
QGP cylinder with a radius Ra- They may dissociate into uncorrelated QQ pairs due to 
multiple scatterings. In this way, the production of heavy quarkionia is suppressed. In the 
following, we will calculate the rapidity dependence of the nuclear modification factor Raa 
in this simplified model. 
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FIG. 5: Heavy-quarkonia suppression by final-state multiple scatterings in most central Au + Au 
collisions: the background medium is assumed to be thermalized within a short time ro ~ 1 fm 28, 



29] , which is still larger than tb , the formation time of heavy quarkonia. After the thermalization 
time To, a quarkonium created at the impact parameter b travels in the QGP cylinder with a radius 
Ra- It may dissociate into uncorrelated QQ pairs due to multiple scatterings if I, the length of 
QGP through which it travels, is long enough. 

A. The nuclear geometry 

In most central collisions, a QGP cylinder with a radius Ra is formed after thermalization. 
We will take Ra to be approximately equal to the radius of the colliding nuclei. Let us 
introduce the nuclear thickness function T^fe) 



T A (b) = J dzp A (z,b), (30) 
where the Woods-Saxon nuclear density pa is given by 

/ \ Pnm fni\ 

PA{r) = — r , w v 31 

1 + exp [(r - R A ) /a\ 



with p nm = 0.17 fm 3 , Ra = 6.38 fm and a = 0.53 fm for Au 24]. The mean number of 
binary collisions in most central collision is defined as 

N coll = J d 2 ba pp Tl(b), (32) 



with the inelastic total cross-section for p + p collisions a pp ~ 40 mb |25l| . 

As illustrated in Fig. [51 after created in the binary nucleon-nucleon collisions at the 
impact parameter b, a heavy quarkonium moves at an angle 9 relative to the cylinder axis. 
Let us calculate I, the length of the background medium through which it travels, i.e., 
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the dashed line protion of the quarkonium's trajectory in Fig. [5j It also depends on the 
azimuthal angle <p. Given (b, 9, <p), one can easily get 

l(b, 6, <p) = ytfA-v^y-bc ostp 

sm 

with l(b, 0, ip) = at b > R A - 



B. The nuclear geometry 

The nuclear modification factor Raa of the heavy quarkonium (QQ) is defined as 

dN^/dy 

RAA = N^dN^/dy' (34) 

with dNQg/dy being the quarkonium (QQ) yield in A + A collisions, N co u the mean number 
of binary collisions in the centrality bin, and dN^/dy the quarkonium (QQ) yield in p + p 
collisions. If heavy ion collisions could be taken as uncorrelated nucleon-nucleon binary 



collisions of the two colliding nuclei, there would be no heavy-quarkonia suppression, i.e., 
Raa = 1- Both the initial cold nuclear effects, such as gluon saturation effect and final- 
, Ijj] can make Raa different from 1. In this paper, we will focus on final-state 



state effects [7 



multiple scattering as discussed in the previous section. With the following observations 

• (a) the quarkonium (QQ) yield in p + p collisions dN^/ dipdy is independent of (p, 

• (b) 7a (6) has rotational symmetry in the impact parameter space, 
we have 

/ d%a pp Tl(b) dcpd 2 N p Q p Q /difdyP sur (r((b, <p, y)), r ) 



Raa ^ Po : 



N coll dN p Q p Q /dy 

J dbba pp T 2 A (b) J d(pP mr (r(l(b,(p,y),r )), (35) 



Ncoll 

where /, the length of the QGP through which the heavy quarkonium travels, is given in 
Equ. fl3"31 and Pq is the contribution from cold nuclear effects. As a first estimate we will 
take P as the average value of Raa due to cold nuclear effects over rapidity y. 

One observation about J/ty production in p + p collisions at ^snn = 200 GeV by the 
PHENIX experiment at RHIC is that the average transverse momentum squared (p^) at 
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forward rapidity is not significant lower than that at midrapidity|26j. As a first estimate, 
we take the typical transverse momentum of J/\I/ produced in binary p + p collisions to 
be independent of rapidity y. In this case, one can easily get their typical momentum p and 
typical velocity vo 



and 



with 



and 



Pt 

V=—^ 36 
sin 9 



My) [mfa + tf)) \m%Jrfl+(p 2 T ) ] ' (37) 

sm9=- -, (38) 

1 + e 2 y v ' 



{p 2 T ) ~ 3.75 GeV 2 . (39) 
Accordingly, the effective length r(l(b,9,<p)) is given by 



r(Z(M,y)) = -I 2 cos 2 9 



P 2 



1 R 2 A -b 2 sin 2 ip - b cos <pj fl + 



m 



(40) 



and 



r( 7 r B ) = ^t 2 - z 2 (t) = A /t 2 - v 2 t 2 cos 2 

m 



'l - v 2 cos 2 0(t o + 7 tb) = r + [ 1 + ) T B , (41) 

where z(t) is the trajectory of the quarkonium. To avoid extending our formula to the 
non-perturbative regime, we need to define a critical invariant time r c , such that, 

T(r c ) = To = T c , (42) 

with T c the critical temperature for the deconfinement phase transition. And we use a 
modified formula for Raa as follows 

p r r2~x 

Raa ^ -r^- / dbba pp Tl(b) / dipP sur (min {r(Z(6, <p, y)), r c } , r ) . (43) 

^ coll J JO 
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FIG. 6: Fits about heavy-quarkonia suppression at RHIC: because of boost invariance of the 
background medium, Raa is independent of y. Our task is reduced to fitting the data about 
J '/iff Raa vs y in most central Au + Au collisions in Ref. [6] to a straight line: Raa = 0.286 with 



X 



9.4. 



TABLE II: Fit results. 



Po 






T c [GeV] 


TO 


[fm] 


T [GeV] 


q [GeV 2 /fm] 


e [GeV/fm 3 ] 


1.0 


0.5 


0.5 


0.19 


0.5 


- 1.5 


0.299 - 0.247 


0.93 - 0.53 


5.44 - 2.55 




1.0 


0.5 


0.19 


0.5 


- 1.5 


0.301 - 0.248 


0.95 - 0.53 


5.62 - 2.60 


0.4 


0.5 


0.5 


0.19 


0.5 


- 1.5 


0.240 - 0.212 


0.48 - 0.33 


2.25 - 1.39 




1.0 


0.5 


0.19 


0.5 


- 1.5 


0.240 - 0.212 


0.48 - 0.33 


2.28 - 1.40 



In the following, we are interested in the nuclear modification factor for J /iff production 
Raa only in the rapidity range —2.0 < y < 2.0. In this case ,we get the Lorentz factor for 
a J /if) with the typical momentum 7 < 1.88. And this justifies that the formation time 
•yr B < t ~ 1 fm. 
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FIG. 7: The degeneracy of the paramters go and To: because Raa is independent of y, we can not 
break the degeneracy of on and To by only fitting the data about J '/^f Raa vs y in most central Au 
+ Au collisions in Ref. [6J. Given to, -Po and £, we can decipher the thermodynamical properties 
at to by fitting RHIC's data as shown in Fig. [6l One can see that the transportation coefficient % 
dose not sensitively depends on our choice of tq in the range 0.5 fm < tq < 1.5 fm and £. 

C. Discussion of the fits 



We use Equ. (j43p and fl25|) for our fits and perform numerically the z integration in 
Equ. (J3DD and the and 6 integration in Equ. (Hgj) . According to (HDj) . (Hill and (|4"31) . 

is independent of rapidity and this is the manifest of Bjorken boost invariance of the 
background medium. Finally, our task is reduced to fitting the data about J Raa vs y in 
most central Au + Au collisions in Ref. jfi] to a straight line. As showed in Fig. [6j in most 
central collisions the best fit for Raa versus y is Raa = 0.286 with x 2 — 9-4- 

Our fit results are presented in Table [IT] and Fig. [71 In Fig. [7J we take q as a function 
of To- To see the effects of the cold nuclear matter to our estimates, we show two sets of 
results with P = 1.0 and Po — 0.40 respectively. In the case with P Q = 1.0 we neglect the 
contribution due to cold nuclear effects. Therefore, we can get the largest estimated value 
of go in our estimate. In the other case with Pq = 0.4 we estimate the cold nuclear effects 
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by taking (pj) = 0? with Q 2 S = Cp/N c 1.25 GeV 2 , the quark saturation momentum squared 
at RHIC's energy 27]. This gives the lower bound to the value of go i n our analysis. The 
exact value of Pq at RHIC's energy requires more careful analysis and beyonds the scope 
of our paper. With £ and P$ fixed, for each value of rn[28L |29[ , we can get a unique value 
for go by fitting RHIC's data. One can see that the transportation coefficient go dose not 
sensitively depends on our choice of To in the range 0.5 fm < To < 1.5 fm and £. In Table 
HI] we list the value of each parameter in our formula. We choose T c ~ 190 MeV, which is 



indicated by the lattice simulation 



301 ] . Our results of g — 0.33 — 0.95 GeV 2 /fm are in sharp 



contrast with those obtained by the analysis of high-p^ hadron spectra in Ref. 3l| in which 



the authors conclude that the values of the time-averaged q exceed 5 GeV 2 /fm. In Ref. 32], 
go is estimated to be 1.6 — 2.0 GeV 2 /fm in central Au+Au collisions at RHIC, which is still 
at least 2 times larger than our results. Therefore, the puzzle that the J/\I/ suppression is 
observed to be smaller than expected from the analysis of the energy density at RHIC as 
discussed in Ref. jsf also shows up in our calculation. 

At the end of this subsection, let us discuss briefly the validity of eikonal approximation 
within At ~ tf orm in describing J/ty production at RHIC energy. Our calculation in Sec. 
Ill CI is justified only if g 2p c ~ pj/y where p is the typical momentum of the (anti- 
)heavy quark and g is the typical momentum transfered under each individual scattering 
between the (anti-)quark and one constituent of the QGP. In a thermal medium, one has 



q~T < T ~ 0.25 GeV. At RHIC energy, 2p c ~ p J/qi > ^/[pj) ~ 1.936 GeV. As a result, 
2^- < 0.13. Therefore, at RHIC energy, the eikonal approximation should be good enough 
for us to estimate the significance of the final-state multiple scattering to heavy-quarkonia 
suppression within At ~ tf orm . 



IV. CONCLUSION 



Due to their small size, heavy quarkonia act as an excellent hard probe to the properties of 
the background QCD matter in heavy ion collisions. In this paper, we focus on their another 
static property, the short formation time tb- If the background medium is thermalized at 
time tq > t^, one can use heavy quarkonia to probe the thermodynamic properties at time 
To- To illustrate our point, we present a simplified model in Fig. [5] to discuss the dissociation 
of heavy quarkonia due to final-state multiple scatterings in Bjorken's expanding QGP. From 
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FIG. 8: The attraction between the quark and anti-quark: and are respectively the relative 
transverse coordinates of the quark and antiquark in the amplitude and in the conjugate amplitude. 
The quark and anti-quark interact with each other via thermal gluons with a mass 

fit results of J/^f Raa versus y in most central Au + Au collisions at A /s/v r /v = 200 GeV at 
RHIC, we get the transportation coefficient go — (0.33 — 0.95) GeV 2 /fm (accordingly the 
energy density eo — (1.39 — 5.62) GeV/fm 3 in perturbative thermal QCD). Our results about 
the transportation coefficient go are in sharp contrast with those obtained by the analysis of 



high-px hadron spectra in Ref. 31]. 



In this work, we use Bjorken's boost-invariant hydrodynamical model to describe the 
evolution of the background medium and obtain a rapidity-independent Raa- Obviously 
this is one reason for the deviation of the measured J/ty Raa at RHIC from such a rapidity- 
independent behavior just as the deviation from the central-plateau structure of particle 
multiplicity in the rapidity spectra observed at RHIC. One refinement of this work is to 
introduce the centrality dependence of Raa by using more complicated hydrodynamic model 
to describe the background medium. Another refinement of this work is to investigate in 
details the competition between the multiple scattering effect and the attraction between the 
quark and antiquark in a quarkonium within a time interval of At ~ Ra ^> Tb- The inclusion 
of the diagrams shown in Fig. [8] in Equ. ( Tl9|) enables one to compare the multiple scattering 
effect, Debye screening and the melting if the imaginary part of fi 2 D included 33j. Besides, 



other possibly important effects have not been addressed, including the gluon saturation 



effects[9j and other final-state effects, e.g., jet quenching |l 31]. the nuclear absorption and 
recombination effects 0, Ijj]. Moreover, in this paper we only discuss the lowest cc bound 
states. In a more quantitative analysis, one should also include the contribution from excited 
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quarkonium states to ground-state production. Accordingly, still another refinement of this 
work is to include all other effects mentioned above to make a quantitative comparison. 

V. ACKNOWLEDGEMENTS 

We are greatly indebted to Prof. A. H. Mueller for numerous illuminating discussions 
and comments. We would also like to thank Prof. R. Baier, Prof. M. Laine, Dr. C. 
Marquet, Dr. A. Vuorinen and Prof. Y. Mao for useful discussions. This work is supported 
by the National Natural Science Foundation of China (Nos. 10721063, 10975003). BW. 
also acknowledges partial financial support from the Humboldt foundation through its Sofja 
Kovalevskaja program. 



[1] K. Adcox et al. [PHENIX Collaboration], Nucl. Phys. A 757, 184 (2005) 

|arXiv:nucl-ex/0410003l . 
[2] An updated review see, R. Rapp and H. van Hees. larXiv:0 903.1096 [hep-ph]. 
[3] T. Matsui and H. Satz, Phys. Lett. B 178, 416 (1986). 
[4] C. Baglin et al. [NA38 Collaboration], Phys. Lett. B 220, 471 (1989). 



[5] H. Satz, arXiv:hep-ph/9806319 



[6] A. Adare et al. [PHENIX Collaboration], Phys. Rev. Lett. 98, 232301 (2007) 
[arXiv:nucl-ex/0611020] . 



[7] R. L. Thews and M. L. Mangano, Phys. Rev. C 73, 014904 (2006) |arXiv:nucl-th/ 0505055 1 . 
[8] L. Grandchamp, S. Lumpkins, D. Sun, H. van Hees and R. Rapp, Phys. Rev. C 73, 064906 

(2006) [arXiv:hep-ph/0507314] . 
[9] D. Kharzeev, E. Levin, M. Nardi and K. Tuchin, Nucl. Phys. A 826, 230 (2009) 
|arXiv:0809.2933l [hep-ph]]; Phys. Rev. Lett. 102, 152301 (2009) [ arXiv: 0808 .2954 [hep-ph]]. 
[10] X. Zhao and R. Rapp, Eur. Phys. J. C 62, 109 (2009) |arXiv:0810.4566l [nucl-th]]. 
[11] Y. Liu, Z. Qu, N. Xu and P. Zhuang, arXiv:0907.2723l [nucl-th]. 
[12] D. Das [STAR Collaboration], Eur. Phys. J. C 62, 95 (2009). 
[13] U. A. Wiedemann. [arXTv:0908.2306l [hep-ph]. 



20 



[14] A. Capella, L. Bravina, E. G. Ferreiro, A. B. Kaidalov, K. Tywoniuk and E. Zabrodin, Eur. 

Phys. J. C 58, 437 (2008) [arXiv:0712.433T1 [hep-ph]]. 
[15] J. D. Bjorken, Phys. Rev. D 27, 140 (1983). 

[16] B. B. Back et al, Phys. Rev. Lett. 91, 052303 (2003) [arXiv:nucl-ex/0210015| . 
[17] R. Baier, Y. L. Dokshitzer, A. H. Mueller and D. Schiff, Phys. Rev. C 58, 1706 (1998) 
[arXiv:hep-ph/9803473] . 

[18] R. Baier, Y. L. Dokshitzer, A. H. Mueller, S. Peigne and D. Schiff, Nucl. Phys. B 484, 265 

(1997) [arXiv:hep-ph/9608322] . 
[19] F. Dominguez and B. Wu, Nucl. Phys. A 818, 246 (2009) |arXiv:0811.1058l [hep-ph]]. 
[20] F. Dominguez, C. Marquet and B. Wu, Nucl. Phys. A 823, 99 (2009) |arXiv:0812.3878l [nucl- 

th]]. 

[21] H. Satz, J. Phys. G 32, R25 (2006) |[arXiv:hep-ph/0512217] . 

[22] R. Baier and D. Schiff, JHEP 0609, 059 (2006) |arXiv:hep-ph/0605183] . 

[23] M. L. Miller, K. Reygers, S. J. Sanders and P. Steinberg, Ann. Rev. Nucl. Part. Sci. 57, 205 

(2007) [arXiv:nucl-ex/0701025| . 
[24] B. Hahn, D. G. Ravenhall and R. Hofstadter, Phys. Rev. 101, 1131 (1956). 
[25] C. Amsler et al. [Particle Data Group], Phys. Lett. B 667, 1 (2008). 

[26] A. Adare et al. [PHENIX Collaboration], Phys. Rev. Lett. 98, 232002 (2007) 

[arXiv:hep-ex/0611020] . 
[27] A. H. Mueller, Nucl. Phys. A 715, 20 (2003) [arXiv:hep-ph/0208278l . 

[28] P. F. Kolb, U. W. Heinz, P. Huovinen, K. J. Eskola and K. Tuominen, Nucl. Phys. A 696, 197 
(2001) [arXiv:hep-ph/0103234] , P. Huovinen, P. F. Kolb, U. W. Heinz, P. V. Ruuskanen and 
S. A. Voloshin, Phys. Lett. B 503, 58 (2001) | arXiv:hep-ph/0101136l , P. F. Kolb, P. Huovinen, 
U. W. Heinz and H. Heiselberg, Phys. Lett. B 500, 232 (2001) |arXiv:hep-ph/0012137| . 

[29] D. Teaney, J. Lauret and E. V. Shuryak, Phys. Rev. Lett. 86, 4783 (2001) 
[arXiv:nucl-th /001 10"58] . 

[30] P. Petreczky, Nucl. Phys. A 830, 11C (2009) [arXiv: 0908 .19171 [hep-ph]]. 

[31] K. J. Eskola, H. Honkanen, C. A. Salgado and U. A. Wiedemann, Nucl. Phys. A 747, 511 
(2005) [arXiv:hep-ph/0406319] . 

[32] F. Arleo, JHEP 0609, 015 (2006) [arXiv:hep-ph/0601075] . 

[33] M. Laine, Nucl. Phys. A 820 (2009) 25C |arXiv:0810.iTT2l [hep-ph]]. 



